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Tracking  Performance  of  the  RLS  Algorithm 
Applied  to  an  Antenna  Array  in  a  Realistic 
Fading  Environment 

Brian  C.  Banister  and  James  R.  Zeidler,  Fellow,  IEEE 


Abstract — In  this  paper,  frequency  domain  techniques  arc  used 
to  derive  the  tracking  properties  of  the  recursive  least  squares 
(RLS)  algorithm  applied  to  an  adaptive  antenna  array  in  a  mobile 
fading  environment,  expanding  the  use  of  such  frequency  domain 
approaches  for  nonstationary  RLS  tracking  to  the  interference 
canceling  problem  that  characterizes  the  use  of  antenna  arrays  in 
mobile  wireless  communications.  The  analysis  focuses  on  the  effect 
of  the  exponential  weighting  of  the  correlation  estimation  Tdtcr 
and  its  effect  on  the  estimations  of  the  time  variant  autocorrelation 
matrix  and  cross-correlation  vector.  Specincally,  the  ease  of  a 
flat  Rayleigh  fading  desired  signal  applied  to  an  array  in  the 
presence  of  static  interferers  is  considered  with  an  AR2  fading 
process  approximating  the  Jakes’  fading  model.  The  result  is  a 
mean  square  error  (MSE)  performance  metric  parameterized  by 
the  fading  bandwidth  and  the  RLS  exponential  weighting  factor, 
allowing  optimal  parameter  selection.  The  analytic  results  arc 
verified  and  demonstrated  with  a  simulation  example. 

Index  Terms — Antenna  array,  fading,  multipath,  recursive  least 
squares,  RLS,  spectral  analysis. 


1.  Introduction 

The  use  of  mobile  wireless  systems  has  been  growing 
rapidly,  and  with  this  growth  has  come  the  need  for 
greater  network  eapaeity  and  a  eorresponding  increase  in  the 
sophistieation  of  the  eommunieation  systems  themselves.  An 
important  area  for  the  enhaneement  of  network  capacity  is 
the  use  of  antenna  arrays,  whieh  ean  provide  performanee 
enhaneement  through  both  interferenee  suppression  and  spatial 
fading  diversity  [1],  [2].  The  mobile  wireless  ehannel  provides 
speeial  difficulties  for  reeeive  array  algorithm  design,  giving 
rise  to  fast  fading  ehannels  that  must  be  aeeurately  estimated 
and  traeked  for  eoherent  modulation  systems.  This  paper 
analyzes  the  tracking  performance  of  such  an  antenna  array 
when  reeursive  least  squares  (RLS)  is  used  to  generate  the 
eombining  weights. 
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Due  to  their  relative  simplieity,  algorithms  employing 
weighted  linear  eombining  of  the  signals  from  eaeh  antenna 
are  most  commonly  considered,  and  many  such  algorithms 
exist  [3].  Wiener  combining  techniques,  which  were  originally 
applied  to  the  antenna  problem  by  Widrow  et  al.  in  [4],  provide 
the  linear  output  estimate  of  the  desired  signal  with  minimum 
mean  square  error  (MMSE).  Wiener  combining  can  be  consid¬ 
ered  optimal  in  the  additional  sense  of  providing  the  maximum 
a  posteriori  (MAP)  signal  estimate  in  the  presence  of  Gaussian 
noise  and  interference  [5],  which  is  commonly  applied  to 
the  decoding  of  received  signals  [2].  Previous  results  on  the 
performance  of  Wiener  algorithms  applied  to  antenna  arrays 
have  derived  the  average  signal-to-interference  ratio  [6]  and 
outage  probabilities  [7]  for  array  systems  with  perfect  tracking 
in  a  fading  environment.  Metrics  such  as  mean  squared  error 
(MSE)  and  bit  error  rate  (BER)  for  receive  antenna  arrays  in 
a  fading  environment  have  been  found  by  simulation  [8]-[10]. 
This  paper  goes  further  in  that  it  provides  analytic  results  for  the 
tracking  performance  of  RLS  in  such  an  environment,  allowing 
optimization  of  performance  in  the  Rayleigh  fading  channel. 

RLS  [11],  typically  exponentially  weighted  (EW-RLS),  is  a 
commonly  considered  adaptive  Wiener  technique.  Tracking  of 
adaptive  prediction  with  RLS  has  been  analyzed  for  zero  band¬ 
width  [12]  and  finite  bandwidth  [13]  chirped  signals.  Optimiza¬ 
tion  of  the  windowing  function  for  prediction  of  a  chirped  signal 
is  considered  in  [  1 4].  Tracking  in  system  identification  problems 
has  been  studied  in  [15]-[18]  for  a  first-order  Markov  process. 
These  previous  results  consider  the  performance  of  RLS  for  pre¬ 
diction  and  identification  problems  in  nonstationary  environ¬ 
ments  with  simple  models  defined  by  a  single  parameter  and  a 
static  or  deterministically  time- variant  receive  vector  autocorre¬ 
lation  matrix.  These  models  do  not  apply  to  the  interference  can¬ 
celing  problem  characterizing  a  receive  antenna  array  in  com¬ 
munications,  and  such  closed-form  tracking  analyses  cannot  be 
found  for  the  interference  canceling  problem. 

Spectral  techniques  have  been  applied  to  RLS  tracking  anal¬ 
ysis  of  system  identification  applications  in  [19]  and  [20].  An 
important  advance  to  the  analysis  of  RLS  tracking  in  [  1 9]  is  the 
derivation  of  the  MSE  for  a  general  class  of  signals  in  terms  of 
the  power  spectral  density  (PSD)  of  the  nonstationarity.  In  addi¬ 
tion  to  deriving  the  MSE  performance  of  EW-RLS,  [19]  shows 
that  weighting  functions  other  than  the  exponential  weighting 
give  better  performance  and  derives  the  optimal  weighting  as  a 
function  of  the  PSDs  of  the  nonstationary  process  and  the  noise. 

This  paper  derives  the  MSE  at  the  output  of  an  RLS  combiner 
for  an  antenna  array  operating  with  a  model  of  a  dynamic  fading 
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Fig.  1 .  Block  diagram  of  ihc  adaptive  combiner  with  crossover  detector  modeling  bit  decisions  for  BER  simulations. 


channel  appropriate  to  mobile  wireless  communications  prob¬ 
lems.  The  system  bloek  diagram  is  shown  in  Fig.  1 .  Similar  to 
[19],  spectral  analysis  tools  are  applied  to  the  analysis  problem. 
However,  it  is  found  that  the  array  problem  introduces  an  addi¬ 
tional  hurdle  to  the  analysis.  In  the  system  presented  in  [19],  the 
autocorrelation  matrix  of  the  received  vector  is  time  invariant, 
whereas  in  the  array  problem,  this  matrix  is  time  variant,  and  the 
resulting  effeets  must  be  carefully  eonsidered.  A  new  technique 
of  quasideterministie  approximation  of  the  least  square  (LS) 
sample  eorrelation  matrix  that  eaptures  its  relationship  with  the 
sample  cross-eorrelation  veetor  is  introduced.  It  is  found  that 
the  effects  of  the  LS  weighting  function  and  fading  PSD  can  be 
separated  from  the  statistical  expectations,  providing  a  tractable 
analytic  solution.  The  results  allow  for  the  selection  of  the  op¬ 
timal  exponential  ‘‘forgetting”  factor  for  EW-RLS. 

In  the  interest  of  tractability,  the  channel  model  has  been  sim¬ 
plified  to  a  level  that  allows  demonstration  of  the  relevant  per¬ 
formance  measures.  Signals  are  represented  in  discrete  time. 
All  impinging  signals  are  considered  to  be  planewaves  from  the 
far-field,  and  the  modulation  is  considered  to  be  “narrowband” 
in  the  sense  that  the  time  delay  across  the  array  can  be  mod¬ 
eled  as  a  carrier  phase  shift.  The  desired  signal  is  degraded  by 
flat  Rayleigh  fading  with  an  approximate  Jakes  Doppler  spec¬ 
trum  defining  the  fading  process,  whereas  the  interfering  signal 
channels  are  static.  While  the  system  under  study  is  assumed 
to  be  DS-CDMA,  issues  of  code  phase  tracking  and  resolv¬ 
able  multipath  are  ignored  in  favor  of  clearly  demonstrating  the 
channel-tracking  properties  of  the  RLS  algorithm.  It  is  worth 
noting,  however,  that  the  RLS  algorithm  would,  in  practice,  have 
to  be  applied  independently  to  each  of  several  resolvable  mul¬ 
tipath  signals  in  a  DS-CDMA  system.  Thus,  the  analysis  here 
presented  would  apply  equally  to  each  such  path. 

II.  Preliminaries 

A.  Vector  Signal  Models  for  a  Signal  and  Interference 

An  overview  of  models  that  have  been  considered  for  use  with 
antenna  arrays  is  presented  in  [21  ].  This  study  uses  the  simplest 
possible  fading  model  and  neglects  issues  of  time  and  angular 
dispersion  in  order  to  focus  on  the  performance  of  the  adaptive 
tracking.  The  ehannel  model  eorresponds  to  that  shown  in  Fig.  2, 
where  the  desired  mobile  experiences  local  multipath  reflec¬ 
tions.  The  angular  dispersion  y?  of  the  desired  mobile’s  signal 
as  reeeived  by  the  base  station  is  small;  therefore,  the  fades  ex¬ 
perienced  by  the  multiple  antennas  arc  identical.  This  channel 


Fig.  2.  Environmental  conditions  modeled  in  the  analysis. 


model  allows  for  performance  enhancement  from  the  array  due 
to  beam  steering  and  null  forming  but  not  from  receive  diver¬ 
sity.  Similarly,  the  time  dispersion  is  small  so  that  multipath  is 
not  resolvable  in  the  receiver.  If  the  local  scatters  are  uniformly 
distributed  around  the  mobile,  then  the  fading  envelope  has  a 
correlation  defined  by  the  Jakes’  model  [22]. 

There  are  M  -f  1  users,  with  baseband  signals  represented 
by  the  (A/  -f  1)  x  1  vector  x[n].  These  modulation  streams  arc 
assumed  to  have  constant  unit  power,  to  be  independent  from 
eaeh  other,  and  to  be  uncorrclatcd  from  time  sample  to  time 
sample  so  that 

£(x[n]x^^[n  —  in]  =  I6[m])  (1) 

where  ^[r?i]  is  the  Kroncckcr  delta  function.  This  signal  model 
may  represent  DS-CDMA  sampled  at  the  chip  rate  or  a  narrow 
band  system  sampled  at  the  symbol  rate. 

a:o[n]  is  the  signal  of  the  desired  transmitter.  It  is  convenient 
at  several  points  to  consider  the  aggregated  interferers,  and  a 
tilde  is  used  throughout  for  this  purpose  so  that 

x[n]  =  [n]  X2[n\  ...  X;i/[n]]^.  (2) 

The  signal  reeeived  by  the  N  antennas  is  a  linear  eombination 
of  the  signals  from  the  M  +  1  users,  as  determined  by  the  time- 
varying  channel  matrix  C[n],  plus  a  white  noise  term  v[n]. 

u[n]  =  C[n]x[n]  -f  v[n].  (3) 

The  channel  matrix  is  comprised  of  the  response  of  the  desired 
signal  and  of  the  interferers,  giving 

C[n]  =  [coN  C]. 


(4) 
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The  white  noise  term  has  an  autocorrelation  given  by 

E{v[n]v^[n  —  m])  =  all6[m].  (5) 

For  an  array  operating  on  signals  from  point  sources  in  the  far 
field,  the  channel  matrix  can  be  given  as 

C[n]  =  DA[n]  (6) 


where  D  is  the  matrix  comprised  of  the  uniform  linear  array 
(ULA)  steering  vectors  [2 1  ]  of  the  M  + 1  users  and  N  antennas, 
and  A[n]  is  the  time-varying  diagonal  matrix  giving  the  time- 
varying  complex  channel  gains  for  each  user.  We  then  have  (7), 
shown  at  the  bottom  of  the  page,  and 


■«oN 


A[71]  = 


«1 


=  SdevoHr) 


OM  J 


(8) 

(9) 


where  ao[n]  is  the  complex  channel  gain  of  the  Rayleigh  fading 
desired  signal,  and  for  z  0  is  a  constant  channel  gain  asso¬ 
ciated  with  each  of  the  interferers.  The  use  of  constant  channel 
gains  for  each  interferer  contributes  to  the  tractability  of  the 
analysis  since  the  value  of  7,  which  was  introduced  in  (78),  is 
constant.  For  a  number  of  interferers  significantly  larger  than 
the  number  of  antennas,  one  would  expect  that  the  results  with 
constant  gains  would  approximate  the  result  with  fading  inter¬ 
ferers.  A  is  the  wavelength  normalized  linear  antenna  spacing 
for  a  carrier  frequency  a;,  and  Oi  is  the  angle  of  arrival  of  the  zth 
signal. 

(jj 

A  = - ; - r  •  (antenna  spacing  distance.) 

27r.  (speed  of  light)  ^  ^  ^  ^ 

For  simplicity  of  notation,  Ec{  •)<=>£(•  |  C[n])  will  denote 
expectation  given  the  random  channel.  The  received  vector  has 
a  time-varying  instantaneous  autocorrelation  matrix  given  by 


R[n]  =  £c(u[n]u^^[n]) 

=  R  +  co[n]c^[n]  (10) 


where  the  eorrelation  matrix  of  the  noise  and  interferenee  is 
given  by 

R  =  CC"+a^I.  (II) 


The  instantaneous  cross-correlation  of  the  received  vector  and 
the  desired  signal  is 


p[n]  =  £c(u[n]a;S[n]) 

=  co[n].  (12) 


B.  Time-  Varying  Channel 

It  will  be  important  to  be  able  to  decompose  the  channel  gain 
of  the  fading  signal  into  terms  that  are  dependent  and  indepen¬ 
dent  of  the  instantaneous  realization  of  the  channel  at  a  given 
point  in  time.  This  decomposition  is  possible  for  a  Gaussian 
random  process  as  the  decomposition  can  be  done  along  the 
lines  of  correlated  and  uncorrelated  components.  This  can  be 
applied  to  the  fading  channel  gain  as  the  Rayleigh  random 
process  (the  gain  magnitude)  arises  from  a  two-dimensional 
(2-D)  Gaussian  complex  gain. 

Let  p[m]  be  the  normalized  autocorrelation  function  of  the 
channel  gain;  then 

£;(ao*[n]ao[n  +  m])  =  aQp[m]  (13) 

In  particular,  p[m]  may  be  the  Bessel  function  of  the  first  kind  of 
zero  order  that  arises  from  Jakes’  model.  Note  that  the  normal¬ 
ization  by  the  mean  fade  power  gives  p[0]  =  1.  The  normalized 
power  spectral  density  (PSD)  of  the  fading  process  is 
which  is  the  Fourier  transform  of  p[n\. 

The  stochastic  channel  gain  can  be  decomposed  into  terms 
relative  to  the  realization  at  time  n,  giving 

ao[7i  +  m]  =  p[rM]ao[n]  +  u,^[n  -F  m].  (14) 

The  complex  Gaussian  process  i>'n[k]  is  defined  relative  to  ao[n], 
and 


E{a[n\j/n  *  [A:])  =  0  (15) 

^O^ri[A:])  =  0  (16) 

*[A:])  =  ~  k]  -  p[i  -  n]p  *[A:  -  n]).  (17) 


III.  Analysis 

A,  RLS  System  Definition  and  Analysis  Approximations 

The  excess  MSE  due  to  weight  misadjustment  is  decomposed 
into  contributions  from  weight  lag  and  weight  noise  [11].  The 
lag  component  arises  from  the  response  of  the  weights  to  the 
time-varying  channel,  whereas  the  noise  component  arises  due 
to  noisy  estimation.  The  general  LS  formulation  is 


oo 

R[n]  =  /z[n  —  (18) 

k=—oo 

oo 

p[n]  =  ^  /i[n  -  ^']u[A:]xo*  [A:]  (19) 

k=—oo 

w[n]  =  R“^[7i]p[n]  (20) 

a:o[n]  =  w^^[n]u[n].  (21) 


1 

j27rA  co5(^o) 
cos(^o) 


1 

^-j27rAcos(^l) 

^-j47rAcos(^i) 


1 

^-j2nA  cos(^ay) 
g— j47r.A  cos(^M  ) 


g-y(2A' -2)7rA  cos(^o)  g-i(2.V-2)7rA  cos(^i  ) 


g-j(2A'-2)7rAcos(^AD 


E>  —  [do  •  •  •  d./\/]  = 


(7) 
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The  sequence  h[n]  is  herein  considered  an  estimation  filter 
that  sums  to  unity  (unit  gain  at  zero  frequency)  and  has  a  Fourier 
transform  given  by  The  EW-RLS  procedure  with  an 

exponential  forgetting  factor  A  <  1  recursively  generates  the 
weighted  sample  autocorrelation  matrix  R  and  cross-correla¬ 
tion  vector  p  with  adaptation  rate  defined  as  =  1  —  A,  where 
larger  ft  gives  a  faster  adaptation.  As  in  [19],  the  performance 
using  the  a  priori  weights  is  considered.  Defining  the  a  priori  es¬ 
timation  filter  impulse  response  h[n\  for  normalized'  EW-RLS 
as 


h[n]  = 


0, 

(1  -  A)A”-^ 


the  a  priori  estimation  error  is 


n  <  0 
n  >  0 


e[ii]  =  ioN  -  a:oN- 


(22) 


(23) 


The  instantaneous  MSE  given  the  channel  and  the  estimates  is^ 

|e[n]|^R[n],p[n]} 

=  p^^[n]R“^[n]R[n]R“^[n]p[n] 

-  p^HR“^[n]p[n] 

-  p^^  [n]R“^  NpN  +  1  (24) 


where  R[n|  is  the  instantaneous  statistical  autocorrelation  ma¬ 
trix,  which  is  still  a  random  function  of  the  random  channel.  In 
order  to  proceed  with  the  tracking  analysis,  the  sample  autocor¬ 
relation  matrix  will  be  assumed  to  reasonably  approximated  by 
considering  it  to  be  quasideterministic,  given  the  instantaneous 
channel 


Rolnj  =  E^(R[n]  |  c[n  —  D\) 

^  n[n]  (25) 

R-^[n]R[n]R-^[n]  ^  E{n[n]  \  c[n  -  D]y^ 

N  (26) 

where  D  represents  some  delay  relative  to  the  realization  time 
n  introduced  by  the  causal  estimation.  This  assumption  is  an 
extension  from  similar  assumptions  that  have  been  frequently 
made  in  previous  tracking  analyses  of  RLS  systems  [16], 
[1 1]-[13],  where  it  has  been  assumed  that  the  sample  corre¬ 
lation  matrix  estimate  is  exactly  the  instantaneous  statistical 
autocorrelation  matrix.  The  approximation  used  herein  is  more 
precise  in  that  effects  of  the  estimation  filter  h[ii]  and  the 
channel  correlation  can  be  considered. 

The  cross-correlation  estimate  can  be  considered  to  be  a  per¬ 
turbation  from  the  true  value  of  the  instantaneous  cross-correla¬ 
tion. 


p[n]  =  Co[n]  +  f[n]  (27) 


'The  summations  arc  normalized  by  ( 1  —  A)  to  conform  to  the  unit  zero  fre¬ 
quency  gain  condition,  emphasizing  the  estimation  performed  and  simplifying 
portraying  the  error  vector/matrix.  The  two  factors  cancel  in  (20);  therefore, 
there  is  no  net  effect. 

^Notc  that  .E(j’o[r/]j'3[7i])  =  1. 


Applying  the  quasideterministic  approximations  and  taking  the 
expectation  over  the  estimates,  the  resulting  expression  is^ 

Jc[n]  =  l-c^^[n]R£,  [n]co[n]+tr  ^R^,  [n]£:c(fNf^^H) j 

(28) 

where  tr( )  is  the  trace  function,  which  is  the  sum  of  the  diagonal 
components  of  a  matrix. 


B.  Sample  Correlation  Matrix  Approximation  Rd[^i] 

The  sample  autocorrelation  matrix  approximation  (25)  cap¬ 
tures  the  effect  of  the  filtering  of  the  “squared”  outer  vector 
product.  Since  each  sample  vector  u[n]  contributes  only  non- 
negatively^  to  any  eigenvalues  of  the  resultant  matrix  R[n],  this 
filtered  version  is  less  likely  to  attain  nulls  than  the  instanta¬ 
neous  realization  of  R[n]  itself,  which  is  important  when  the 
fading  channel  gain  term  a[n]  in  the  inverted  matrix  may  be 
canceling  the  channel  gain  term  in  p[n].  Further,  the  estima¬ 
tion  filter  h[n\  has  a  nonzero  filter  delay,  and  one  would  expect 
that  the  estimate  of  R[n]  generated  will  thus  have  a  certain  lag. 
The  approximation  defined  here  allows  for  some  consideration 
of  these  effects.  It  is  clear  that  the  expectation  of  the  interferer 
terms  of  (25)  can  extracted  from  the  summation,  leaving 

R£?[ri]  =  R  -h  I  h[n  —  A:]c[A:]c^^  [A:] 

\A=^oo 

(29) 

The  channel  vector  is  then  decomposed  into  components  that 
are  correlated  and  uncorrelated  with  the  fading  process  at  time 
(n  —  £)),  according  to  the  decomposition  of  (14),  as  is  done 
elsewhere  in  the  lag  analysis  of  Section  Ill-E. 


co[n  -  D] 


RdM  =  R  +  dd"  •  £  ^  /|[71  -  fc]|p[^-  -n  +  D] 


fc=— oo 


X  a[n  -  D]-\-  iyn-D[f^]\^ 


a[n  —  D]  .  (30) 


Taking  this  expectation  leaves 

Koin]  =  R -h  .  (  \a[n  -  D]\^  ^  h[n  -  k] 


h=~-oo 

oo 


\p[k  -  n D]\^ -\- Oq  ^  li[n  -  k] 


k=- 


X  (l-\p[k-n^D]\^)]. 


(31) 


Define  the  summation 

oo 

h[-kMk  +  D]\^.  (32) 

h=—oo 

^Hcrc,  the  matrix  relations  tr(AB)  =  tr(BA)  and  a"Ba  =  tr(a"Ba)are 
used,  as  used  repeatedly  throughout  the  analysis. 

^Note  that  this  statement  requires  that  the  filter  response  h[n]  ehosen  is  non¬ 
negative  as,  for  example,  in  standard  EW-RLS  or  rcetangular  windowed  LS. 
This  is  required  for  strietly  ensured  non-negative  definiteness  of  the  resultant 
estimate. 
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Equivalently,  this  can  be  portrayed  in  the  frequency  domain  as 


where  0  denotes  cyclic  convolution.  Note  that  the  integral  sub¬ 
scripts  of  /x,y  are  selected  throughout  to  refer  to  the  order  of 
the  filter  response  x  and  the  PSD  y,  respectively,  in  (33)  and 
(56)-{58).  The  delay  argument  D  appears  only  when  the  order 
of  h  is  unity  (i.e.,  and  Ii^i)  so  that  there  is  no  temporal 
symmetry. 

Then 

RdH  =  R  +  dd''  •  (\a[n  -  +4{^-  4^^))  • 

(34) 

The  inverse  of  this  matrix  is  considered  in  Appendix  A. 

The  purpose  of  incorporating  the  delay  component  D  is  to 
portray  the  sample  correlation  matrix  R[n]  in  a  simplified  form 
that  captures  the  relationship  between  this  matrix  and  the  sample 
cross-correlation  vector  p[n],  and  the  value  of  D  is  selected  with 
this  consideration.  Setting  D  =  0  \s  the  most  simple  choice 
but  ignores  the  delay  effects  of  the  estimation  filter;  another 
simple  but  inadequate  selection  might  be  the  dominant  filter 
group  delay,  which  for  EW-RLS  is  given  by  A/(l  —  A)  samples 
at  zero  frequency.  In  order  to  portray  R[n]  as  function  of  a  single 
instantiation  of  a[n  —  D]  and  best  capture  the  relationship  be¬ 
tween  R[n]  and  p[n],  the  delay  should  be  chosen  to  maximize 
the  residual  relationship  after  taking  expectation.  This  is  best 
accomplished  by  selecting  D  to  maximize  the  cross-correlation 
between  the  delayed  a[n  —  D]  of  the  correlation  matrix  and  the 
filtered  a[n]  of  the  cross-correlation  vector,  which  is  given  by 
[defined  in  (57)].  Hence,  D  is  selected  to  maximize  l[^i  . 

It  is  notable  that  this  approximation  to  the  time-varying 
sample  autocorrelation  matrix  reduces  to  the  instantaneous 
statistical  autocorrelation  when  D  =  0  and  Ii^2  —  ^  » 

the  estimation  filter  is  replaced  with  a  delta  function  in  time), 
providing  an  approximation  similar  to  that  of  previous  works 
[11]-[13],  [16];  however,  as  we  will  discuss  later,  this  ap¬ 
proximation  would  be  inadequate  to  accurately  describe  the 
performance  of  the  present  system. 


C.  Characterizing  RLS  MSE  in  Terms  of  Lag  and  Noise 

Let  the  error  involved  in  each  sample  cross  correlation  be 
z[n],  which  consists  of  the  noise  and  interference  terms.  It 


idefincd  as  follows: 


z[n]  =  xo*[7i]  •  +  v[n]j  .  (35) 

The  first  and  second  moments  of  the  cross  correlation  sample 
error  are  then  given  by 

Ec{z[n\)  =  0  (36) 

Ec{z[n]z^[n  -  r/i])  =  R(5[m].  (37) 

The  lag  and  noise  terms  of  the  sample  cross  correlation  vector 
error  arc  defined  as 

oo 

PiagH  =  'Ej  “  ^']co[^']  -  Co[n]  (38) 

A:=— OO 

oo 

Pnoise[n]  =  h^7i  A:]z[A:]  (39) 

k-=—oo 

and  the  error  vector  is  given  by 

f^n]  —  Plag[^]  P noise  [nj*  (^^) 


The  lag  and  noise  error  vectors  arc  independent  so  that  the  au¬ 
tocorrelation  of  f  is 


£c(fNf"[n])  =  £^c(PlasNPla«"N) 

4" noise [77] P noise  (41) 

In  order  to  obtain  an  overall  cost  measure,  the  expectation  is 
performed  over  the  channel  variation.  This  final  cost,  averaged 
over  the  channel  variations,  is 


J  =  E{Jc[n]) 


—  «7opt  “h  *^IaK 


lag  ■ 


(42) 


From  (28),  (38),  and  (39),  the  components  of  J  arc  found  to 
be  as  in  (43)-(45),  shown  at  the  bottom  of  the  page.  It  is  of 
interest  to  compare  these  results  with  those  of  [19].  As  in  [19], 
the  noise  and  lag  components  of  the  cross-correlation  error  can 
be  considered  to  arise  from  a  noise  filter  /i noise M  (which  is 
referred  to  as  /i « in  [  1 9])  and  a  lag  filter  /lUg  [n]  (which  is  referred 
to  as  hte  in  [19]),  and  the  summations  are  convolutions  with 
these  filter  responses.  These  two  filters  produce  the  combining 
weight  errors  from  noise  and  adaptation  delay  (“lag”).  From 
(38)  and  (39),  these  impulse  responses  are  given  by 

^Lioisefri’j  —  /l[7?J  (^6) 

/iiag[n]  =  /i[7i]  -  6[n].  (47) 


^pt  =  1  -£(co''[n]R  ^[7i]co[7i]) 
j|ag  =  £(c(/^[7i]R“^[n]co[7i])  +  £  ^  tr  I  R^  [n]  ^  ^  /t[n  -  ^-]/7  [n  -  7]co[7]co^^[fc] 


OO  OO 


t  =  — oo  A:  =— oo 


i  -1 


J  -i 


-£|tr|^Ro  [ti]  ^  /i[7i -7]co[7]c5'[n]j  I  -  £|tr  ^Ro  [tx]  ^  /7[n-7]co[7i]c^^[7] 

‘^noise  =  £  I  ir  (  Rd  N  ~  *]  I  I  • 


i=— oo  A:=— oo 


(43) 

(44) 

(45) 
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Thus,  with  0  here  denoting  time  convolution 

f[n]  =  ftiagH  0  Co[n]  +  ftnoiseh]  0  zN-  (48) 

In  contrast  to  the  system  identification  analysis  of  [19],  this 
result  is  not  yet  definitive  since  in  the  analysis  of  [19],  the  au¬ 
tocorrelation  matrix  R[n]  is  explicitly  assumed  to  be  time  in¬ 
variant  based  on  the  statistics  of  the  source.  In  the  context  of 
system  identification,  this  is  a  reasonable  assumption.  However, 
this  condition  is  clearly  not  valid  for  the  interference  cancella¬ 
tion  application  now  under  analysis.  The  MSE  cannot  be  char¬ 
acterized  using  only  these  filter  values  because  there  is  an  inter¬ 
action  between  the  time  variation  of  the  correlation  matrix  and 
the  cross-correlation  vector.  This  is  made  explicit  when  consid¬ 
ering  that  Rd^]  and  R[n]  in  (43 are  both  functions  of 
the  sequence  Co[n],  which  also  appears  in  the  sample  cross-cor¬ 
relation  error  summations.  Thus,  the  present  problem  requires 
a  more  rigorous  consideration  of  the  interplay  between  these 
functions.  This  is  presented  in  the  following  sections. 


D.  Optimal  Time  Varying  Wiener  Combiner 

The  MSE  performance  of  the  time  varying  Wiener  combiner 
with  no  tracking  error  was  given  in  (43).  Using  the  matrix  in¬ 
version  lemma  result  (77)  from  Appendix  A  gives 

l«oN|2 


^opt 


l-E 


7  +  l<afoH|" 


(49) 


where  ao[n]  is  complex  Gaussian  so  that  |cvoMP  follows  an 
exponential  distribution  with  mean  ctq,  and 
112  .  .  _  . 

+  1 


and 


4  I 

V7  +  I«oNIV 

^pt  =  -^e^/"»Ei(-7/a2) 


(50) 

(51) 


where  Ei( )  is  the  exponential  integral  function  [23]. 


E.  Excess  MSE  Due  to  Lag 

The  lag  term  (44)  of  the  MSE  is  the  most  complicated  as 
it  incorporates  the  interaction  of  the  channel  variation  in  the 
cross-correlation  vector  and  in  the  autocorrelation  matrix. 

Recall  that  the  angular  dispersion  of  the  arriving  multipath  is 
assumed  to  be  very  small  (see  Fig.  2),  giving  rise  to  correlated 
fading  across  all  of  the  antenna  elements.  In  this  case,  there  is 
no  diversity  gain  (only  beam  forming  gain).  This  simplification 
allows  us  to  decompose  the  correlation  matrix  and  apply  the  ma¬ 
trix  inversion  lemma  as  in  Appendix  A.  Then,  using  the  results 


of  (77)  with  (34)  but  separating  cvo[*],  where  it  cannot  be  ex¬ 
tracted  from  the  summation,  (44)  becomes 


'^lag  —  ^ 


(  l«oMP  ^ 

U+l«o[n]|V 


f  T,Z-oo  l>b>-  -  -  A:]ao[z]ao*[A;l\ 

V  7  +  <^o  (l  -  4?)  +  |(»[n  -  j 


I  T,Z-oo  -  ^]Re(ao[^]ao•  [n])  \ 


Then,  cto[i]  is  decomposed  relative  to  time  (n  —  D),  as  in 
(14),  such  that  the  expectation  can  be  separated  along  the  lines 
of  independent  components,  as  in  (53),  shown  at  the  bottom  of 
the  page.  Extracting  ao[n  —  D]  from  the  summation,  taking  the 
expectation  within  the  summation  using  the  autocorrelation  of 
the  orthogonal  process  i^n[  ']  (17),  and  discarding  the  vanishing 
cross  terms  gives  (54),  shown  at  the  bottom  of  the  next  page. 
The  form  of  (54)  shows  that  the  effects  of  the  channel  gain  PDF 
have  been  decoupled  from  the  effects  of  the  channel  gain  auto¬ 
correlation.  Further  analysis  is  then  broken  into  finding  the  ex¬ 
pectation,  which  characterizes  the  channel  PDF,  and  computing 
the  summations,  which  contain  the  channel  correlation  function. 
A  similar  decoupling  could  be  found  for  Ricean  or  Nakagami 
fading  PDF  as  they  can  both  be  characterized  by  Gaussian  pro¬ 
cesses. 

Of  the  expectations  in  (54),  one  was  already  considered  (50), 
and  the  other  is  given  by 


E 


(D) 


1  f  1  + -  A', 

- exp  - - i - 


'*1,2  ^0 


AD)  2 
'*1,2  ^0 


X  Ei  - 


7  +  ag  (l 


•*1,2  ^0 


Define  inverse  transform  integrals  for  7i  =  0 


’  27r  J_^ 


(55) 


(56) 

(57) 


,  _p/  Mn]|^ 

Y,Z-oo  ~  ^  +  -Dlaoh  -  £>]  +  ;^„-D[j])(p*[^'  -  n  +  -Dloo'h  -  77]  +  ;/„-D*[fc]) 


+ 


t(d) 


7  +  <^o"(l-/{°^)+|a[n-I7]|2/{^ 

2E^-oo^[”  “  ’']Rc((p[i  -n  +  Z7]ao[n  -  Z7]  +  i^„_D[i|)(p*[i7]ao*[n  -  D]  +  )7.-d.M)) 


7  + 


(53) 
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l2,o=^  (58) 

Applying  the  results  for  the  expectations  (50),  (55),  and  the 
frequency  domain  integrals  {56)-(58)  yields 


F.  Excess  MSE  Due  to  Weight  Noise  Error 
The  excess  MSE  due  to  weight  noise  error  was  given  in  (45). 

Again,  decomposing  Il[n  —  D]  and  applying  the  matrix  inver¬ 
sion  lemma  gives 


J.,g  =  +  1 

^0 


Equation  (59)  provides  the  MSE  contribution  due  to  lag,  incor¬ 
porating  the  channel  fading  correlation  function  in  the  integrals 
(56)-(58)  so  that  the  lag  MSE  can  be  quickly  obtained  for  dif¬ 
ferent  mobility  rates  or  channel  assumptions. 


CO  oo 

'^noise  ~  ^  ^  ^  ^  ll^Tl  ~~  Aj]5[z  —  Aj] 

i=— CO  k=—oo 

_ |ao[n  -  D]\'^ _ 

7  +  ag  (l  -  4°^)  +  |a[n  - 

CO 

X  ^  ^  h[n  —  z]/i[7i  —  k]6[i  —  A:],  (60) 

i=— oo  k=—oo 


Through  similar  procedures  to  those  used  for  the  lag  calcula¬ 
tion  and  using  (50) 


l2,0  T'  +  *^0  (l  - 


'^noise  -  (o) 

■‘1,2 


X  exp 


X  Ei  - 


tW  2 
-^1,2  ^0 

7  +  ag  (l  - 


;(D)  2 
'*1,2  ^0 


(61) 


From  (6 1 ),  it  is  seen  that  the  PSD  of  the  fading  channel  docs  not 
affect  the  noise  misadjustment  term,  although  the  PDF  of  the 
channel  gain  is  incorporated. 


-Alai?  =  E 


l«oNP 

7  +  |o[n]|2 


+  E 


\ao[n-D]\‘^ 


CO  CO 


y  7  +  (l  -  +  l«[n  -  J  J 

X  p[z  —  n  +  D]p*  —  n  4-  D] 

1 


-  iVi[n  -  k] 


CO  OO 


+  aoE 


-alE 


+ <^0  (i  -  4a)  + 1«[«  -  w4a  I  ooii-f: 


h[n  —  i]h[n  —  k]p[i  —  A;] 


oo  CO 


V7  +  (l  -  4^^)  +  |o[n  -  j  ii 

X  /j[i  —  n  +  D\p'  [^•  —  n  +  D\ 

A  \ao[n-DW _ 


-2E 


+  (i  -  ^[^2)  + 


(D) 

2  , 


Re  I  p*[D]  •  ^  h[n  —  i]p[i  -  n  +  D] 


2oqE 


+  2<7q£’ 


'  +  (1  - 1[^4)  +  l«h  -  *= 


/i[n  —  z]Re(p[n  —  ?]) 


h[n  —  i]Re(/i*  [D\p[i  —  n  +  D]). 


(54) 


1044 


IEEE  TRANSACTIONS  ON  SIGNAL  PROCESSING,  VOL.  50,  NO.  5.  MAY  2002 


Fig.  3.  Comparison  of  results  for  Jakes  model  (69)  and  the  AR2  model  (71)-(74)  with  {D  —  0,  fixed  /[®2  =  with  (D  =  Dopt,  =  .r,  denoting 

normal  ealeulation).  ealeulated  numerically  for  Jakes’  spectrum.  Doppler  frequency  of  fa  =  Ir'  —  4  cyclcs/samplc.  Shown  for  varying  SNR. 


G.  Summary  of  MSE  in  Rayleigh  Fading 

The  MSE  of  the  RLS  antenna  array  system  has  been  decom¬ 
posed  into  the  Wiener  MSE,  excess  MSE  due  to  lag  misadjust- 
ment,  and  excess  MSE  due  to  noise  misadjustment.  The  overall 
MSE  of  the  system  is  given  by  the  summation  of  these  three 
MSE  terms,  which  is  given  by  (5 1 ),  (59),  and  (6 1 ).  The  analysis 
has  allowed  for  the  separation  of  the  complex  Gaussian  fading 
gain  PDF  from  the  PSD  of  the  fading  process.  The  effeets  of  the 
time  correlation  of  the  fading  process  arc  captured  in  the  four 
integrals  (33),  (56),  (57),  and  (58).  The  “optimaP’  value  of  D 
is  termed  Z^opt  and  selected  to  maximize  as  discussed  in 
Section  I1I~B.  The  MSE  is  calculated  as 


—  Z)opt  — 


J  —  1  -{- 


=  arg  (max  |) 

-  2Re  -  hfi 


AD) 


S  +  ag  (l  -  /l',2  j 


- 


AD)  2 
^1,2  ^0 


exp 


'7  + (To  (1  - 


'*1,2  ^0 


X  Ei  -- 


7  + 


(1 


+ 


r(D)| 

^1,1 


-h 


2Re  -  P*[D]l[^f  -  /2j 


r(.D) 

'1,2 


X  exp 


AD)  2 

^1,2  ^0 


X  Ei  - 


7  + 


,2 ; 


n 


(D>„2 


artificially  fixed  =  1.  The  Jakes’  PSD  is 

2 

PjakesCe-^'")  =  ' 


Cj|  <  CJd 


(64) 


0,  |cj|  >  u}d 

/^jakes[^l']  =  'Zo(u^<i7l)  (65) 

where  u)d  is  the  Doppler  frequency,  and  Jo{x)  is  the  zero-order 
Bessel  function  of  the  first  kind. 

For  the  filter  of  (22)  and  PSD  of  (64),  the  integrals  (56)  and 
(57)  can  be  found  with  Z)  =  0  for  small  by  using  the  approx¬ 
imation  z  ^  1  ju,  and  where  appropriate,  A  ^  1. 

ff  2 


^2,1 


(62) 


r 

ft 


r  du) 


/(o)  r 
27r  J_, 


IP  -  jpu 


w,i  (P 


c/cj 


=  ^2,1. 

The  final  integral  (58)  is  given  by 

(l-A) 


(66) 


(67) 


72,0  = 


(1  +  A) 
2-ir 


(68) 


Hence,  with  fixed  /}*2  =  H  the  MSE  of  the  system  is  given  by 
7|d=o  =  (-7/<t;^) 

V^o 


-h 


- 1 


+  1 


(63) 


+ 


ft 


fft 


s/lP  + 


ZP,  (P  +  w, 


(-et^^'oEi  (-7^0)) 


ft 


■  (  ^c‘'7‘'oEi  (-7/a^ 


IV.  Numerical  Examples 
A.  Jakes  ’  Model  Fading 

For  the  Jakes’  model  [22],  approximate  closed-form  results 
are  available  using  the  simplified  case  where  D  =  Q  and  with  eral  case  Z)  ^  0  and/or  realistic  be  done  numerically. 


0  —  R  V  iT^2  \  '  /  ■-'u  /  /  •  (69) 

2  P  \^o  J 

This  MSE  result  is  compared  with  that  of  the  AR2  model  (see 
Section  IV-B)  in  Fig.  3.  Evaluation  of  the  integrals  for  the  gen- 

(^), 
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Fig.  4.  Autocorrelation  of  Jakes’  model  and  AR2  approximation  /d  =  O.OOOl  cyclcs/sample. 


B.  AR2  Fading 

In  order  to  provide  tractable  mathematical  expressions  for  this 
numerical  example,  the  Jakes’  model  will  be  approximated  as  a 
seeond-order  autoregressive  process  (AR2),  with  poles  selected 
to  provide  a  PSD  close  to  that  of  the  Jakes’  model.  This  model 
provides  not  only  an  analytic  solution  but  a  simple  implemen¬ 
tation  for  simulation.  A  good  match  between  the  Jakes’  correla¬ 
tion  and  the  AR2  correlation  for  small  time  offsets  is  found  by 
selecting  the  poles  according  to  the  rule 

j)  =  ^  j  exp(±j  •  0.8  •  LJd)  (70) 

where  ujd  is  the  discrete  time  Doppler  frequency,  and  the  poles 
are  in  the  ^-plane.  The  AR2  process  and  the  resulting  PSD  is 
described  in  more  detail  in  Appendix  C. 

By  way  of  a  specific  example,  consider  a  system  with  a  1 
MHz  signal  bandwidth  and  100  Hz  fading,  as  would  arise  from 
a  velocity  of  67  mi/h  (108  km/h)  with  a  carrier  frequency  of  1 
GHz.  In  the  discrete  time  domain  after  sampling  at  the  1  MHz 


Nyquist  rate,  these  parameters  give  a  fading  rate  of  le  —  4  cy¬ 
cles/sample  or  cjt/  =  27r  •  le  —  4  rad/sample.  The  AR  model  is 
then  defined  with  poles  at  0.9998  •  exp(±j27r  •  0.000  08).  The 
resultant  autocorrelation  is  shown  in  Fig.  4.  Simulations  are  also 
run  with  /</  —  5c  —  4  cyclcs/sample  for  comparison. 

For  an  AR2  process  with  poles  at p  and/?*,  the  integral  values 
are  given  in  (71)-(74),  shown  at  the  bottom  of  the  page,  using 
the  AR2  nomenclature  of  Appendix  C. 

Plugging  (71)-(74)  into  (63)  provides  a  closed-form  value  of 
the  MSE  at  the  output  of  the  antenna  array  combiner. 

C.  Simulation  and  Discussion 

Simulations  were  performed  with  the  AR2  model  in  order 
to  confirm  the  analysis.  The  simulation  results  are  compared 
to  several  analytic  configurations,  which  show  that  the  delay 
term  D  must  be  properly  selected  as  described  in  (62),  and 
cannot  be  simplified  to  unity  but  rather  must  be  explicitly  deter¬ 
mined.  Since  ultimately  a  metric  such  as  bit  error  rate  is  of  in- 


(D)  ^  (l-A)AP  /A-^r^-A-^-^r^^+^  1 

l-h  005(2-0)  \  1  —  1  — 

A'"^*^  cos(20  -F  26)  —  ^-^-17*2^5+2  cog^20  +  2D6  -(-  26)  —  cos(20)  +  cos(2D^  -h  20) 

'''  1  -  2A-ir2  cos(26l)  + 


•^2,1  — 


I 


_ 

1,1  “ 


>^2,0  — 


co5(20)  —  Ar^  cos(— 2^  +  20)  \ 

^  1  —  2Ar2  cos(2^)  -h  A^r**^  / 

_ (l-A)-(l^r^)-|l-p^|sin(^) _ 

(1  -F  A)(l  -pA)(A  -  p){l  -p*A)(A  -p*)sin(^  -  angle(l  -p^)) 
(1  -  A)^  sm(2(9  -  angle((l  -  p^){\  -  Ap)(A  -  p)))  2 

|1  —  Ap||A  — /?|  sin(^  —  angle(l  —  p^)) 

'  sin(-D(9  +  26  +  angle((l  -  Ap)(l  -  p^))) 

<  (l--^)s^nW(l-^)|l~r^l 

co.s(v)(A^  — 2Arco«(^)4- T'^)(l  — 2X7*  cos(d)+A^7-^) 

+  ^  -  angle((A  -  p)(l  -  p^))) 
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Fig.  5.  MSE  from  analysis  (solid  lines)  and  simulation  (points).  Doppler  frequeneies  —  4,  uc  —  4  eyeles/saniple,  as  indieated. 
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Fig.  6.  Cross-over  probability  (BER)  from  simulation.  Doppler  frequeneies  le  —  4,5r— 4.  Perfeet  traeking  lines  are  aeeording  to  Gaussian  interferenee 
approximation  with  fading  eonstant  throughout  a  bit  symbol. 


terest,  BER  was  extracted  from  the  simulation  in  addition  to  the 
MSE.  The  simulation  implemented  a  pseudo-CDMA  system,  as 
shown  in  Fig.  1 ,  where  the  signal  is  known  by  both  receiver  and 
transmitter  in  order  to  facilitate  tracking,  as  is  the  case  when  a 
training  sequence  or  pilot  is  available.  Since  the  signal  is  known 
in  this  case,  there  are  no  information  bits  to  extract.  To  extract  a 
BER  metric,  bit  decisions  are  modeled  as  a  crossover  decision. 
For  CDMA  reception,  the  signal  is  multiplied  by  a  code  (x*[n]) 
and  accumulated.  In  this  case,  the  “data"  is  all  +ls;  therefore,  if 
the  accumulator  output  is  greater  than  zero,  then  there  is  no  de¬ 
tection  error,  and  if  it  is  less  than  zero,  then  there  is  a  detection 
error.  This  can  equivalently  be  considered  to  be  a  system  with 
perfect  decision-directed  feedback.  This  provides  an  adequate 
model  for  a  general  evaluation  of  the  tracking  performance  of 
the  system. 

It  is  assumed  that  each  transmitter  is  power  controlled  so  that 
the  average  power  received  by  the  base  station  from  each  mo¬ 
bile  is  the  same,  but  the  desired  signal  mean  power  (Tq  is  varied 
over  [1  10  100]  to  better  illustrate  the  BER  performance.  The 
interferers  experience  a  static  channel  with  unit  mean  received 


power,  and  the  desired  signal  undergoes  Rayleigh  fading.  The 
antenna  array  is  comprised  of  three  elements,  spaced  0.4  wave¬ 
lengths  apart.  The  desired  transmitter  has  an  angle-of-arrival 
(AOA)  of  108°,  and  there  are  eight  interferers  with  randomly 
assigned  AOAs  of  32°,  73°,  111°,  133°,  143°,  157°,  165°,  166°, 
and  1 68°.  No  AWGN  is  applied  to  the  received  vector;  therefore, 
all  noise  considered  is  due  to  the  interferers;  this  emphasizes  the 
gains  from  null  forming. 

With  these  parameters,  the  effective  signal  to  noise  is  0.125  • 
<7q  with  one  antenna,  whereas  with  the  three  antennas,  one  can 
obtain 

£(|ao[n]|2)  •  d^'R-Mo  =  O.oUal 

This  is  a  6.1  dB  improvement  over  the  single  antenna  case. 

The  simulations  were  run  using  AR2  fading  rates  of  Ic— 4  and 
5c— 4  cycles/samplc.  The  results  are  seen  in  Figs.  5  and  6,  which 
show  the  mean  square  error  and  bit  error  rate  of  the  system  as  a 
function  of  the  adaptation  rate  ft.  Noting  that  ft  is  the  3  dB  band¬ 
width  in  units  of  rad/sample  of  the  estimation  filter,  the  x-axes 
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Fig.  7.  Group  delay  of  cslimation  filter  H )  with  ,.i  =  —  3  (3  dB  bandwidth  (2;r)  ^  —  3  eyelcs/sample). 


Fig.  8.  Time-varying  combiner  weight,  illustrating  the  weight  delay,  =  l,^o/7  =  —2.9  dB,  =  Ir  —  iji  =  \r  —  3,  delay  =  999  samples.  Real 
component  of  iri  is  shown. 


arc  scaled  by  l/27r.  The  MSE  simulations  match  the  analysis 
for  smaller  /i,  becoming  less  accurate  for  larger  /i,  where  the 
noisy  sample  autocorrelation  matrix  degrades  performance  and 
the  quasidcterministic  assumption  breaks  down. 

Both  the  analysis  and  simulations  indicate  that  in  order  to  at¬ 
tain  reasonable  performance,  the  3  db  bandwidth  arising  from  ft 
must  be  selected  to  be  about  an  order  of  magnitude  larger  than 
the  Doppler  frequency  (Figs.  5  and  6).  To  consider  why  so  much 
margin  is  needed  in  fi,  consider  the  filter  response  The 

delay  of  the  estimation  filter  h[n]  at  low  frequencies  is  (1  -  ft)/f3 
samples,  shown  in  Fig.  7.  This  delay  is  seen  to  be  very  signif¬ 
icant  near  the  transition  band  when  considered  in  units  of  de¬ 
grees.  Thus,  the  large  margin  in  the  value  of  ft  is  required  to  ac¬ 
commodate  this  delay,  which,  otherwise,  will  cause  lag  error  in 
the  eombincr.  A  plot  of  one  of  the  time-varying  weights,  along 
with  the  optimal  weight  and  a  delayed  version  of  the  optimal 
weight,  is  given  in  Fig.  8.  As  the  problem  formulation  is  similar 
aeross  the  weights,  showing  one  weight  is  adequate  for  quali¬ 
tative  visualization.  This  figure  clearly  shows  that  for  moderate 
lag  errors,  the  lag  takes  the  form  of  a  minimally  distorted  but 
delayed  version  of  the  optimal  weights,  and  the  zero  frequeney 
filter  delay  of  (1  —  which  is  999  samples  in  the  figure, 
eharacterizes  this  delay  well.  Most  of  the  fading  process  power 


is  near  the  edge  of  the  band,  where  the  group  delay  effects  arc 
strongest,  exacerbating  this  problem.  The  plot  of  delay  in  de¬ 
grees  (see  Fig.  7)  shows  that  about  one  decade  of  margin  be¬ 
tween  Doppler  frequency  to  the  filter  3  dB  bandwidth  is  required 
in  order  to  avoid  this  lag  effect.  It  is  also  of  interest  to  note  that 
the  BER  performance  of  the  be  —  A  eyelcs/sample  Doppler  case 
is  actually  lower  for  higher  SNR  than  the  le  -  4  eyelcs/sample 
case  when  /^  =  27r  •  le  —  2,  counter  to  expectations  in  a  tracking 
problem;  this  is  due  to  time  diversity  as  the  very  deep  fades 
that  produce  bit  errors  at  high  SNRs  arc  also  very  brief,  and  at 
faster  fading  rates,  the  fades  do  not  annihilate  the  entire  64-chip 
symbol. 

Fig.  9  shows  MSE  simulated  and  calculated  with  various 
values  of  D  and  l[^\  demonstrating  that  the  formulation  of 
the  complete  analysis  is  required  for  accurate  results.  With  the 
“optimal”  D  and  calculated  the  results  closely  match  the 
simulation.  The  simplest  form  to  calculate  would  be  D  =  0  and 
artificially  set  =  1;  this  form  gives  large  peaking  when 
the  Doppler  frequency  falls  in  the  estimation  filter  transition 
band.  Letting  D  =  Dopt  and  setting  =  1  also  results  in  a 
large  error.  The  error  in  these  cases  is  particularly  large  for  high 
SNRs  because  here,  the  eorrcct  handling  of  the  cancellation 
of  the  a[n]  terms  in  the  inverted  autocorrelation  matrix  R  and 
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Fig.  9.  MSE  calculations  for  various  values  D  and  <7^  =  100.  Dopt  is  the  value  of  D  that  maximizes  D^c  =  A/(l  —  A)  is  the  filter  group  delay 
at  zero  frequency;  =  r  indicates  normal  calculation;  =  1  fixes  the  value  to  unity. 


the  cross-correlation  vector  p  becomes  more  important.  The 
filtering  present  in  the  sample  autocorrelation  matrix  R  means 
that  even  when  the  signal  fades,  there  is  residual  filtered  signal 
in  the  sample  autocorrelation  matrix  so  that  near  divide-by-zero 
conditions  arc  avoided  in  the  realized  R;  this  motivates  the 
definition  of  R,  which  leaves  a  “lagged  residual”  of  a[n]  in  the 
form  of  the  ((1  —  term.  It  is  also  clear  from  the  figure 

that  the  selection  of  the  value  of  D  is  important. 

The  result  for  the  AR2  model  [see  (71)-(74)]  is  compared 
with  the  result  for  the  Jakes’  model  [see  (69)]  in  Fig.  3  with  the 
simplifying  conditions  D  =  0  and  fixed  /|  2  =  discussed 
above,  these  conditions  do  not  allow  a  very  accurate  prediction 
of  the  performance,  but  the  closed-form  approximate  result  for 
the  Jakes’  model  is  available.  Results  obtained  by  numerical 
integration  for  the  Jakes’  model  are  also  displayed  for  the  most 
accurate  formulation  [see  (62)  and  (63)].  The  calculated  MSE 
for  the  Jakes’  and  AR2  models  are  nearly  identical,  illustrating 
that  the  use  of  an  AR2  approximation  to  Jakes’  is  appropriate. 

The  performance  is  limited  by  the  exponential  weighting  of 
the  EW-RLS  estimates.  We  have  surmised  that  this  is  largely 
due  to  the  phase  shift  imposed  by  the  exponential  filter.  Lin  et 
ai  [19]  suggested  optimal  weighting  functions  to  replace  the 
exponential  weighting  for  the  identification  problem.  Results 
similar  to  those  could  be  found  here  from  (63).  However,  the 
significance  of  this  would  be  obscured  by  the  complexity  of  the 
equations  involved. 

It  is  of  interest  to  compare  the  present  tracking  results  with 
those  derived  for  system  identification  of  a  Jakes  fading  channel 
[19]  and  the  first-order  Markov  process  [16],  as  well  as  for 
prediction  of  a  chirped  process  [12].  The  noise  misadjustment 
is  found  to  be  approximately  proportional  to  ft  in  all  of  these 
cases,  which  is  in  agreement  with  (74).  The  lag  misadjustment 
in  [  1 6]  is  found  to  be  inversely  proportional  to  ft,  whereas  in  [  1 9] 
and  [12],  it  is  found  to  be  inversely  proportional  to  ft^.  From 
(71)-(74),  it  is  clear  that  no  such  simplification  is,  in  general, 
possible  for  the  interference  canceling  mobile  wireless  commu¬ 
nications  problem.  The  lag  misadjustment  is  a  function  of  the 


fading  PSD  and,  in  general,  cannot  be  reduced  to  some  simple 
polynomial  of  ft. 


V.  Conclusion 

This  paper  has  presented  a  derivation  of  the  MSE  that  results 
from  the  application  of  the  EW-RLS  algorithm  to  an  antenna 
array  receiving  signals  degraded  by  a  Rayleigh  fading  channel. 
It  is  shown  that  the  time  correlation  of  the  fading  channel  can 
be  extracted  from  the  probability  density  function  of  the  channel 
such  that  the  exponential  weighting  used  for  the  cross-correla¬ 
tion  estimate  could  be  considered  to  be  a  filtering  operation. 
Considering  the  windowing  function  to  be  a  filtering  opera¬ 
tion  allows  the  application  of  well-known  Fourier  analysis  tech¬ 
niques  so  that  the  analysis  can  be  conveniently  applied  to  chan¬ 
nels  with  different  fading  PSDs. 

The  application  of  Fourier  analysis  is  similar  to  [19]  but  is 
applicable  to  the  “interference  canceling”  tracking  problem  that 
characterizes  the  antenna  array  problem.  This  paper  expands  the 
use  of  such  techniques,  illustrating  a  methodology  for  applying 
such  analysis  when  the  correlation  matrix  of  the  received  vector 
is  time  varying,  as  is  frequently  the  case  in  communications 
problems. 


Appendix  A 

Autocorrelation  Matrix  Inverse  Decompositions 


A.  Decomposition  of  the  Inverse  of  the  Statistical  Correlation 
Matrix  i2[n] 


It  is  frequently  necessary  throughout  the  analysis  to  separate 
the  contributions  of  the  desired  signal  and  the  interferers  to 
the  autocorrelation  matrix  inverse.  This  can  be  accomplished 
through  straightforward  application  of  the  matrix  inversion 
lemma. 

The  decomposed  correlation  matrix  is  given  by  (10).  The  ma¬ 
trix  inversion  lemma  applied  to  (10)  gives 


R-^n]  =  R-^  - 

1  +  c^^[n]R“^co[n] 


(75) 
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The  inner  product  of  the  channel  response  with  respect  to  the  in¬ 
verse  of  the  correlation  matrix  can  be  reduced  to  a  simple  form. 


Co^[n]R  ^[n]co[n] 


c^^[n]R  ^co[n] 

1  +c5^[n]R"‘^co[n] 


(76) 


Given  the  narrow  angular  dispersion  of  the  channel  model,  this 
reduces  to 


c^[n]R  ^[n]co[n] 


7  +  |ao[n]|2 


(77) 


where  7  provides  a  sense  of  the  irreducible  noise  power  pro¬ 
jected  into  the  signal  space. 

-  =  d'^R-Mo.  (78) 

7 


B.  Decomposition  of  the  Inverse  of  the  Sample  Correlation 
Matrix  Estimate  Rd[^] 

The  matrix  inversion  lemma  applied  to  (34)  gives 

Al\n]  =  R-^  -  (|a[n  -  +  cr^  (l  - 

R-Mpd^^R-^ 

(79) 

The  inner  product  of  the  channel  response  with  respect  to  the  in¬ 
verse  of  the  correlation  matrix  can  be  reduced  to  a  simple  form. 


c"[n  -  i7]RD  [n]co[n  -  D] 


\a\[n-D]\^ 


{D) 

2 


(80) 


Appendix  B 

It  is  informative  to  consider  the  SNR  generated  by  the  optimal 
Wiener  combiner.  In  the  context  of  a  time-varying  system,  this 
provides  the  time-varying  SNR  that  would  result  at  the  output 
of  a  Wiener  combiner  with  perfect  tracking  capabilities. 

The  Wiener  estimator  of  xo[n]  is  given  by 

rro[n]  =  CoHR“^[n]u[n].  (81) 


be  obtained  by  simply  using  the  common  equations  for  fading 
[2].  This  has  been  used  to  provide  the  “perfect  tracking”  bit  error 
rates  of  Fig.  6,  with  the  approximation  of  the  interferer  signals 
as  Gaussian.  It  is  well  known  that  this  approximation  holds  well 
in  DS-CDMA  systems  for  low  SNRs,  but  performance  at  higher 
SNRs  can  be  better  than  indicated  with  this  approximation  be¬ 
cause  the  tails  of  the  binomial  distribution  are  not  as  long  as 
those  of  the  Gaussian  distribution. 


Appendix  C 

Properties  of  the  AR2  Process 

The  AR2  process  of  interest  has  poles  that  form  a  conjugate 
pair  at  p  and  p*,  where  p  =  t*  exp(j0),  0  <  0  <  tt,  excited  by  a 
white  Gaussian  process.  The  generating  filter  under  considera¬ 
tion  will  be  normalized  to  give  an  output  power  unchanged  from 
the  power  of  the  white  input  process. 

The  autocorrelation  function  of  such  an  AR2  process  with 
unit  power  input  can  be  found  in  [11].  With 


tl)  =  e  -  --  angle(l  -  p^) 
D  =  sm{9)  ■  (1  —  7-^)  •  |1  — 


(86) 

(87) 

and  with  the  process  normalized  to  obtain  a  unit  output  power, 
the  resulting  autocorrelation  function  is 

jM 

p[n]  =  - — cos(|7il^  -h  VO  (88) 


and  the  PSD  is 


|c-os(V7)| 
D 


|cos(V-^)| 


1 


(1  -p^-i)(i  -j)*^-i)(i  -p2)(i  -p*zy 
The  power  normalized  two-pole  filter  applied  to  the  white  exci¬ 
tation  is  then 

^  1 


•f^AR2(^)||p|<|z|  = 


D 


(90) 


cos(V0  (i  ~ 

The  squared  autocorrelation  is  required  for  (32).  This  is  given 
by 

7.2|nl 


pM  = 


1  -h  cos(2V-^) 


(l-hcos(2|n|6>-f-2V^)).  (91) 


The  signal  power  at  the  output  (conditioned  on  the  channel)  is 

Ec  (|c"[n]R“^[7i]coHa:o[n]|^^  =  (co^[n]R“^[7i]coM)^  . 

(82) 


The  total  power  of  i’oM  *s 
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